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摘要 : 蒸 散 是 土壤 -植物 -大 气 连 续 体 (SPAC) 水 热 运 移 的 一 个 重要 环节 ,是 全 球 水 量 平 衡 的 重要 组 成 ,一 直 是 气象 学 ,水 文学 .地 
理学 及 生态 学 等 相关 学 科 重 要 的 研究 主题 。 区 域 尺度 地 表 蒸 散 的 时 空 变化 过 程 十 分 复 如 3 迄今 为 止 ;在 像 元 太 度 水 平 上 ,特别 
在 非 均 匀 下 垫 面 和 地 形 起 伏 条 件 下 ,有 代表 性 的 进行 地 表 蒸 散 的 观测 仍然 十 分 困难 ,里 然 遥感 方法 可 获得 区 域 尺 度 水 平 落 散 ， 
但 其 主要 根据 经 验 或 半 经 验 模 型 对 区 域 蒸 散 进行 估算 ,模型 选用 的 参数 以 及 结果 还 需 地 面 实测 数据 进行 改进 、 优 化 ,如 何 获 得 
与 交感 尺度 相应 的 地 面 蒸 散 实测 数据 成 为 模型 验证 的 重点 和 难点 , 光 闪 烁 方法 的 出 现 为 上 述 问 题 的 解决 带 来 了 希望 。 光 闪烁 
方法 能 够 适应 复杂 下 垫 面 , 测 量 结果 准确 且 具 有 时 空 平 均等 优点 ,成 为 测量 区 域 蒸 散 的 有 效 方法 、 验 证 遥感 模型 结果 的 最 佳 手 
段 。 从 理论 原理 .计算 方法 .主要 应 用 情况 等 方面 ,概述 了 光 闪 烁 方法 观测 区 域 蒸 散 研究 进展 指出 了 影响 测算 精度 的 不 确定 性 
因素 ,并 提出 了 研究 展望 , 旨 在 进一步 推动 该 方法 在 区 域 莱 散 观测 研究 中 的 应 用 ,促进 相关 学 科 的 发 展 。 

关键 词 :区 域 平均 葵 散 ;大 孔径 闪烁 仪 ;微波 闪烁 仪 ;遥感 ; 通 量 
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Abstract; Evapotranspiration is key to hydrological and heat transfer in the soil-plant-atmosphere continuum ( SPAC). It is 
also an important. component of water and energy balances, and now it is becoming an important research subject in 
meteorology, hydrology , ecology and other related areas. The temporal and spatial variations in surface evapotranspiration at 
the regional scale are very complicated, which means that it is still difficult to measure surface evapotranspiration at the 
kilometer scale above a heterogeneous underlying terrain. Remote sensing methods can provide regional evapotranspiration 
results by using an empirical or semi-empirical model, but the parameters in the models and the outputs need to be 
improved and optimized according to the real-time measurement results. The difficulty with model verification is how to 
obtain surface evapotranspiration at the kilometer scale. However, the development of the scintillometry method may be able 
to solve this problem. The scintillometry method can adapt to complex surfaces, and produces accurate results when both 


time and space are averaged. This method has become an effective way to measure regional evapotranspiration and may be 
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the best way to verify the outputs of remote sensing models. This study discusses the background theory to the scintillometry 
method , summarizes the potential applications of the method , and analyzes the uncertainties in the hardware equipment, the 
environmental factors, and key parameters in the calculation. Finally, we put forward research prospects based on 
theoretical principle, the calculation method, and potential applications. Our summary and analysis can be used to promote 


the application of the scintillometry method in regional evapotranspiration measurements and relevant disciplines. 


Key Words: regional-averaged evapotranspiration; large aperture scintillometer; microwave scintillometer; remote 


sensing; flux 


2 BUE E HE-IUM AK Ofe BE MOK Tas TE RU CERIS TS ,是 全 球 水 循环 的 重要 组 成 部 分 。 在 全 球 气候 变化 的 
背景 下 ,气温 降水 均 发 生 了 显著 的 变化 ,水 循环 的 加 速 导 致 区 域 水 资源 分 布 不 均匀 的 现象 更 为 严重 5 因此 ， 


程 对 天 气 预 报 .干旱 监测 .区 域 水 资源 管理 及 全 球 变化 等 研究 具有 重要 意义 。 

蒸 散 过 程 复杂 ,涉及 范围 广 , 包 括 中 小 尺度 ( 几 百 米 范围 ) 和 区 域 尺度 (公里 范围 六 其 中 以 波 文 比 -能 量 平 
衡 法 .梯度 迭代 法 空气 动力 学 方法 , 涡 度 相关 方法 以 及 燕 渗 仪 法 等 计算 < 观测 方法 为 主 的 中 小 尺度 燕 散 研究 
成 果 丰 硕 , 近 十 年 来 涡 动 相关 方法 因 其 准确 快速 的 优点 常 被 研究 者 认为 是 蒸 散 观测 的 标准 方法 。 由 于 全 
球 变化 和 流域 水 文 研 究 的 需要 ,区 域 尺度 的 蒸 散 研 究 越 来 越 受 到 重视 “3 然而 公里 尺度 的 区 域 燕 散 过 程 影 
响 机 制 较 多 ”2 , 涡 动 相关 方法 观测 结果 并 不 理想 。 若 将 基于 中 小 尺度 的 单 点 式 测 量 结果 扩展 到 区 域 尺度 ， 
必须 考虑 下 垫 面 的 异 质 性 ,如 植被 的 覆盖 率 地形 的 起 伏 程度 等 ,同时 还 需 考 虑 大 气 状 态 的 平稳 性 汕 流 的 发 
展 以 及 局 地 微 气 象 条 件 的 差异 ,因此 要 获得 区 域 平均 莹 散 需 建立 观测 网 络 !"] 。 近 年 来 ,利用 遥感 信息 估算 大 
尺度 地 表 蒸 散 , 即 遥 感 模型 法 ,成 为 观测 区 域 薰 散 的 有 效 技术 之 一 ,但 该 技术 通常 是 根据 经 验 或 半 经 验 模型 对 
DOG 2E BGIET HET UIT 。 由 于 对 莹 散 过 程 中 许多 关键 要 素 的 参数 化 计算 方法 不 是 十 分 成 熟 ,下 垫 面 的 非 均 
匀 性 导致 下 感 蒸 散 模 型 在 不 同 区 域 的 适用 性 不 同 以 及 遥感 反 演 地 表 参 数 的 不 确定 性 等 ,遥感 模型 结果 必须 结 
合 地 面 观 测 数据 进行 验证 改进 ,模型 参数 也 需要 结合 地 面 实测 数据 进行 优化 2 。 如 何 获得 与 遥感 蒸 散 模 
型 尺度 相应 的 地 面 薰 散 实测 数据 成 为 模型 验证 的 前 提 和 重点 ,也 是 遥感 水 文 监测 的 技术 瓶颈 。20 世纪 70 年 
代 发 展 的 基于 汕 流 大 气 光 传播 理论 与 相似 理论 相 结合 光 闪 烁 方法 为 上 述 问题 的 解决 带 来 了 和 希望 。 

1978 年 Wang 提出 了 利用 光 闪 烁 法 测量 通 量 的 假设 ,美国 NOAA 波 传播 实验 室 将 此 假设 变 成 了 可 
BEUA 。 随 后 的 几 十 年 间 ,基于 近 红 外 光 微波 可见光 的 各 种 光 闪 烁 方法 测量 通 量 的 仪器 迅速 发 展 ,并 于 九 
十 年 代 中 后 期 开始 用 于 实际 观测 研究 。 光 闪烁 方法 主要 分 为 “ 单 波长 方法 ”和 “ 双 波 长 方法 ”。 单 波长 方法 是 
基于 近 红 外 方式 (0:67 一 0:94pm) 的 光 闪 烁 仪器 ,此 类 仪器 主要 包括 小 孔径 闪烁 仪 ,大 孔径 闪烁 仪 和 超大 孔径 
闪烁 仪 ,其 中 -SAS 测量 距离 仅 为 200 一 250m ,XLAS 的 测量 距离 可 达 10km, 目前 光 闪 烁 方法 的 应 用 以 LAS 为 
EU, 与 单 波长 方法 不 同 , 双 波 长 方法 是 通过 红外 -微波 联合 系统 工作 ,该 联合 系统 包括 近 红 外 闪烁 仪 和 微 
波 闪烁 仪 , 目前 .OMS 尚 处 于 初级 阶段 。2000 年 ,中 和 荷 合作 项 目 “ 中 国 能 量 与 水 平衡 监测 系统 ”首次 引进 LAS 
以 来 ,先后 在 黄河 流域 .黑河 流域 海河 流域 .黄土 高 原 . 青 藏 高 原 等 地 方 展开 了 针对 LAS. 的 专门 研究 ,并 成 功 
研制 出 具有 自主 知识 产权 的 闪烁 仪 2225 。 光 闪烁 方法 可 以 测量 区 域 范围 内 的 平均 通 量 , 既 弥补 了 传统 观测 
方法 空间 代表 性 不 足 的 问题 ,同时 也 与 遥感 的 象 元 尺度 匹配 度 高 , 较 好 地 验证 了 迁 感 模型 的 反 演 结果 … 7 。 
本 文 从 光 闪 烁 方法 的 理论 原理 出 发 ,着 重 介绍 该 方法 在 区 域 蒸 散 测 量 领 域 的 计算 方法 ,概述 了 该 方法 测算 区 
域 蒸 散 的 研究 进展 ,分 析 了 该 方法 在 应 用 计算 过 程 中 的 不 确定 性 并 对 其 在 今后 的 研究 热点 进行 展望 , 骨 在 进 
一 步 推动 该 方法 在 区 域 蔡 散 观 测 研 究 中 的 应 用 。 


1 光 闪 烁 方法 测定 显 热 和 潜 热 的 原理 


光 在 消 流 大 气 中 传播 , 由 于 空气 温度 ,湿度 和 压强 的 波动 引起 大 气 折射 率 发 生变 化 ,从 而 导致 光束 的 频 
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AE 振幅 强度 等 出 现 了 无 规则 的 变化 ,产生 了 所 谓 的 “ 光 闪 烁 "现象 "5 。 光 闪烁 方法 就 是 根据 接收 光 强 变 
化 的 信和 号 来 按 公式 (1) 来 测定 空气 折射 指数 结构 参数 ( C?) : 
C as = 1.120, DPL? (1) 
六 过 下 过 万 三 大 (2) 
式 (2) 是 式 (1) 的 限定 条 件 。 式 中 C; us 是 LAS 测算 出 的 空气 折射 指数 结构 参数 (m 77) ; o, 为 接收 光 
强 信 号 的 自然 对 数 方差 ;D 为 闪烁 仪 的 孔径 尺寸 Am ) iL 为 光路 长 度 (m) ;7 为 满 流 结构 内 尺度 (m) i F AIRE 
FER BE (m) iL, gri Viti 3] P JN. BE (m) 。 
C 受 空 气 的 温度 (7) BE (q) UREA CP) 的 波动 而 发 生变 化 ,因此 可 用 温度 结构 参数 ( C7 ) WEBER 
构 参 数 ( C; ) 和 温度 与 湿度 的 协 方差 项 ( Cy ) 表示 如 下 : 
A EL. A? 
= Crt- Cr + Cs (3) 
7 Tq q 
INIP „Ap 4, 分 别 是 温度 .湿度 结构 参数 系数 (/) ; 了 是 平均 温度 (K) ; g 表示 平均 湿度 (kg/kg)。 
在 近 红 外 方式 工作 时 ,温度 的 波动 对 C2 的 影响 最 大 ,47 比 4 大 2 一 3 个 量 级 。 为 了 便于 计算 ,假设 温度 和 
湿度 的 相关 系数 为 1, 并 引入 表示 波 文 比 系数 (6) 对 公式 (3) 进 行 了 人 简化: 


47 0.03 
> 2 à 
C, EO. iex (4) 


2 
C, 


APF, 7 是 空气 温度 (K) ,B 为 波 文 比 系数 (/)。 
在 红外 -微波 双 波 长 工作 方式 下 ,不仅 需要 根据 公式 (1) 计 算 近 红外 波段 的 Ci as ,同时 还 需 根据 公式 (5) 
计算 微波 波段 的 C2 wws : 
C) uws ZOS OÀ (5) 
式 中 , Tes EIC 2/ A 为 波长 ,m) ,为 光路 长 度 (m)y a, 为 接收 光 强 信号 的 自然 对 数 方差 (/) 。 公 式 
(3) 因此 出 现 了 两 种 不 同 波长 形式 表达 。 为 顺利 计算 出 2 , C5 ,假设 C, = + 上 RnwWCrC; , Ry 是 经 验 值 , 常 取 


app 
e — At uiis Caes * Peer, mws + ZR p uwsAr ras C ust ms (6) 
q (IT) 
A uL xe d up A D C. 
p q,MWS 7 n,LAS q,LAS ^n,MWS Tq^^q,MWS^ >q, LAS n,LAS ~ n,MWS 
C= RE (7) 
(Tq) 
A i ATA B EM x er p 
Ap T- EE DORUM. 4 和 4 的 下 标 分 别 表示 在 LAS 与 MWS 工作 波段 下 的 温 湿度 参数 取 值 。 
q 


这 种 根据 假设 计算 的 方法 ,被 称 为 “Hill 双 波 长 法 ”或 “ 双 波 长 假设 法 ”。Lidi SEU 认为 波 文 比 系数 变化 时 ， 
R;, 也 会 发 生变 化 ,假设 Rn =+1 的 理想 条 件 并 非 总 是 成 立 , 因 此 提出 获得 实时 Ri 值 的 交互 方法 , 称 为 “ 双 波 
长 相关 法 ”或 “Liidi 双 波 长 方法 ”。 该 方法 不 仅 可 以 计算 出 C; Las 和 CT mws ,而 且 同 时 计算 出 近 红外 与 微波 交 
互 的 空气 折射 指数 结构 参数 C ous o IP C ons 的 计算 公式 如 下 : 

Chons = Ak "PL Pai, (8) 
式 中 ,4 是 待定 系数 ,与 LAS MWS 的 路 径 长 度 和 其 间距 有 关 , 具 体 计算 方式 可 参考 Lidi, AR) 因此 
可 变 成 3 种 不 同 表现 形式 , 联 立 可 得 Cr 、C; 以 及 C, 如下: 


c ius 
Cy, -M' CA (9) 
C C; ous 
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式 中 , MK” 为 道 和 矩阵 ,参考 表达 式 如 下 : 


A, MWS A, LAS E 2A, 1AsÁ q,MWS 
q 9 g 
Mz Jik g " E: Ar wws Au uws T Ar rasy ns Ar Las Aes ras | Ar wwsA s mws (10) 
(Ar ias A, mws a Ar uwsA, tas ) : Tq Tq Tq 
AT. mws AT Las — 2År LisA v mws 
T T T: 


根据 莫 宁 - 奥 布 霍 夫 相似 理论 (MOST) 存在 如 下 关系 ，; 


Cr(z-d | z-d 
go AL (n) 
C (az - d)? EE 
I aer (12) 
qx MO 


式 中 ,z 为 有 效 高 度 (m) ,d 为 零 平面 位 移 (m) ,7, 为 温度 变量 (K) a 为 湿度 变量 (kg/kg) 1.2 MOST HAR 
数 ,rwo 为 奥 布 霍 夫 长 度 。 其 中 ， 
u, T 
Luo =- mA (13) 
式 中 ,g 为 重力 加 速度 (9.81m/s) u, HERE ( m/s) , K Ivon Karmian 常 数 (0.4)。 由 于 Lyo 受 五 和 ,的 
影响 ,上 述 公 式 通常 无 法 得 出 五 或 LE 的 解析 解 , 对 于 红外 -微波 双 波 长 工作 方式 的 光 办 烁 法 需 结合 下 式 进 行 
迭代 运算 。 


* 


H =- pe u, Py (14) 
EET L (15) 
LE 

E= L (16) 
式 中 ,5 为 空气 密度 (1.2kg/m ) ,6 为 空气 比 热 (1005]A( kg * K) ) ,7 是 绝对 湿度 (g/m ) , LÆRER J 
kg) ,H H WPNE E (W/m ) ,LE 是 潜 热 通 量 (W/m ) , 术 7 是 蒸 散 量 (mm) 。 对 于 近 红 外 方式 工作 的 光 闪 烁 方 

法 只 能 迭代 计算 出 克 , 然 后 根据 能 量 平 衔 方 程 , 利 用 余 项 法 得 出 也 己 : 
LE-R,-G,-H (17) 


式 中 R, 是 净 辐 射 (WAm”);G 是 土壤 热 通 量 (W/m?)。 
MOST 普 适 通 数 通常 分 为 如 下 两 种 : 
不 稳定 状态 ( <0): 


fi) = Xp (1 x d 23 f) — X1 (1 — x, o6) p (18) 
稳定 状态 C& 0) 
ju reu Su UL ur) (19) 
Jh 27 是 稳定 度 函数 ，xr，，xra soa sua 是 MOST 的 函数 系数 。 


此 外 ,也 有 使 用 可 见 光 . 近 红 外 、 微 波 3 个 波段 联合 的 方法 进行 蒸 散 测算 .”” ,虽然 这 些 实验 表明 了 使 用 
3 个 波段 工作 的 可 行 性 ,但 实验 结果 并 不 理想 , 且 技 术 要 求 高 .成 本 昂贵 ,逐渐 失去 了 吸引 力 。 


2 光 闪 烁 方法 观测 蒸 散 的 研究 进展 


21 单 波长 方法 
以 近 红 外 方式 工作 的 单 波 长 光 闪 烁 方法 通常 以 波长 为 880nm 或 940nm 的 LAS 为 主要 观测 仪器 。 该 方法 
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发 展 初期 ,在 平坦 下 垫 面条 件 下 对 其 进行 观测 研究 ,并 与 涡 动 相关 方法 测量 结果 进行 比较 。Purerua 半岛 的 平 
坦 草原 进行 的 多 次 观测 实验 结果 表明 ,LAS 的 测量 结果 与 涡 动 相关 方法 的 结果 较 好 ,标准 误差 在 25W/m 之 
内 ,并 认为 LAS 可 以 在 短 时 间 内 获得 路 径 上 的 平均 通 量 ,其 统计 不 确定 性 较 小 站。Cain 等 ' 引 在 牛津 的 实验 
结果 表明 LAS 与 涡 动 相关 方法 的 观测 结果 一 致 性 较 高 。 同 类 研究 还 包括 Kohsiek 等 (站 利用 XLAS 进行 实验 、 
Hoedjes 等 '55 在 墨西哥 西北 部 实验 以 及 Beyrich 等 :59 在 德国 东北 部 的 实验 等 ,这 些 实验 结果 表明 :LAS 的 测量 
结果 与 涡 动 相关 方法 测量 结果 一 致 性 显著 。 国 内 在 2000 年 首次 引进 LAS 以 来 ,结合 相关 气象 资料 ,分别 对 
RE HEE 郑州 等 地 区 初步 进行 水 热 通 量 的 研究 对 比 ,结果 认为 LAS 测量 的 显 热 通 量 与 当地 的 气温 差 相 关 
性 显著 ,结合 净 辐 射 资料 得 出 的 蒸 散 值 均 能 反映 当地 能 量 的 季节 变化 规律 '”2 。 随 后 黄河 流域 海河 流域 、 
黑河 流域 .黄土 高 原 .青藏 高 原 等 不 同 地 区 的 LAS 应 用 试验 表明 ,LAS 观测 结果 与 涡 动 相关 方法 , 波 文 比方 法 
测量 结果 一 致 性 很 高 ”3 。 其 他 研究 如 在 农田 生态 系统 . 半 干 旱地 区 以 及 干旱 地 区 等 结果 也 表明 光 闪 烁 
方法 测量 区 域 通 量 结果 可 靠 !331。 

实际 研究 中 ,区域 尺 度 的 下 垫 面 大 部 分 是 非 均匀 的 ,由 于 地 形 起 伏 、 植 被 异 质 等 客观 因素 的 存在 , 光 闪 烁 
方法 在 复杂 下 垫 面 条 件 下 的 适用 性 颇 受 关注 。Hemakumara 等 在 斯 里 兰 卡 Horana 的 混合 植被 区 进行 实 
验 ,结果 表明 LAS 所 测 蒸 散 量 符合 当地 的 气候 条 件 ,并 能 随 着 天 气 的 条 件 发 生变 化 ,由 此 拉 开 了 光 闪 烁 方法 
在 非 均匀 下 垫 面条 件 应 用 的 序幕 。Beyrich 等 525] 对 德国 东北 地 区 森林 和 农田 混合 区 域 研究 时 ,将 各 地 段 测 得 
的 通 量 利用 权重 函数 求 和 得 到 区 域 平均 通 量 ,该 结果 与 光 闪 烁 方法 吻合 度 很 高 。- 张 劲松 等 5 在 地 形 起 伏 的 
中 国 北方 低地 山 丘 人 工 混 交 林 的 研究 中 将 两 台 涡 动 相关 方法 的 观测 结果 进行 权重 琶 加 ,所 得 结果 与 光 内 烁 方 
法 测量 结果 的 相关 性 系数 可 达 0.93。Chehbouni 等 的 实验 中 也 得 出 相同 的 结论 2 。 郝 小 浴 等 ' 引 研究 陀 东 黄 
土 高 原 下 热 面 不 均匀 性 对 光 闪 烁 方法 的 影响 ,发 现 地 表 不 均匀 性 对 测量 结果 影响 显著 。 在 非 均 匀 下 热 面 时 ， 
由 于 地 形 条 件 等 因素 的 存在 ,应 用 光 闪 烁 方法 测量 时 ,相似 理论 未 必 成 立 , 但 上 述 非 均匀 下 热 面 的 实验 结果 表 
明 ,考虑 采样 尺度 ,并 利用 权重 函数 对 不 同 地 段 通 量 值 进行 修正 , 光 闪 烁 法 同样 可 以 取得 良好 结果 。 光 闪烁 方 
法 在 非 均 匀 下 垫 面 的 良好 结果 说 明 , 光 闪 烁 方法 在 实际 应 用 中 可 以 与 遥感 模型 进行 对 比 验证 。Watts 5509 48 
据 AVHRR 数据 以 及 LAS 观测 结果 对 墨西哥 北部 半 干 旱地 区 的 地 表 通 量 进行 研究 ,认为 LAS 能 较 好 反映 出 
该 地 区 的 通 量变 化 ;Meijninger 等 5 在 土耳其 西部 灌溉 区 域 用 LAS 分 别 测量 2700m 和 670m 距离 的 通 量 值 ， 
结果 与 LANDSAT 结果 具有 良好 的 一 致 性 ,并 认为 LAS 可 作为 独立 的 区 域 平均 能 量 观测 系统 。Nobuhle 5:09 
对 南非 不 同 生态 气候 区 的 研究 表明 ?闪闪 烁 方法 测量 结果 与 卫星 遥感 模型 运算 结果 一 致 性 高 。 同 样 结 果 也 出 
现在 黑河 流域 青藏 高 原 等 实验 中 路 。 总 体 说 来 ,以 近 红 外 方式 工作 的 光 闪 烁 方法 适应 不 同 条 件 下 热 
面 , 且 该 技术 已 经 成 熟 ; 已 经 实现 了 商业 化 生产 ,其 中 常见 的 有 答 兰 Kipp&Zonen 公司 和 德国 Scintec 公司 。 单 
波长 方法 观测 结果 与 遥感 模型 结果 吻合 度 高 ,已 成 为 检验 卫星 遥感 反 演 模型 的 最 佳 方 法 '"'”1。 
2.2” 双 波长 方法 

红外 -微波 双 波长 方式 工作 的 光 闪 烁 方法 还 处 于 起 步 阶段 。Kohsiek 等 ' 趾 最 先 采 用 波长 为 10mm 发 射频 
率 为 30GHz 的 微波 成 功 计算 出 潜 热 通 量 ,并 认为 C, , C? 与 Cn 存在 某 种 关系 ,甚至 认为 同时 假设 C? 与 
JC; C; fü X DU. C/C 的 取 值 后 ,使 用 单独 波段 即 可 实现 水 汽 通 量 的 计算 。1988 年 Hill 采用 173GHz 的 频 
率 研究 认为 低频 波动 的 水 汽 会 对 结果 产生 影响 ,首次 完整 提出 “ 双 波 长 "方法 计算 过 程 , 并 假设 Cj, = 
£RQACLC, ,其 中 R21 7, MeMillan 等 系统 地 分 析 了 不 同 波长 测量 的 结果 ,并 认为 理想 波长 为 11lmm 
或 3.2mm。 从 对 水 汽 敏 感 程度 以 及 成 本 因素 考虑 ,这 两 种 波长 成 为 以 后 MWS 试验 的 首选 。Liidi 55 9 JA JS 
Wt Ri 的 方法 并 不 能 真实 反映 出 波 文 比 ,并 利用 微波 和 近 红 外 交互 测量 的 方法 得 出 与 涡 动 相关 方法 具有 高 吻 
合 度 的 结果 。Meijninger 等 '2 验证 了 双 波 长 方法 在 非 均匀 下 垫 面 测算 区 域 蒸 散 的 可 行 性 。Ward 4% 
国 Swindon 城郊 区 域 利用 双 波 长 方法 研究 水 热 通 量 ,测量 结果 变化 趋势 明显 ,与 涡 动 相关 方法 结果 一 致 性 较 
高 。 国 内 关于 双 波 长 方法 的 研究 较 少 ,舒婷 "5 对 双 波 长 方法 测量 区 域 蒸 散 进行 研究 ,研制 了 一 套 监 测 范围 在 
600m 以 内 的 设备 ,并 在 中 国 北方 草原 下 垫 面 条 件 下 进行 观测 ,结果 表明 双 波 长 方法 得 到 的 显 热 通 量 与 涡 动 相 
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关 方 法 测量 结果 相关 系数 为 0.96 , 洪 热 通 量 的 相关 系数 达到 0.76 ,表明 该 系统 在 草地 下 垫 面条 件 下 具有 一 定 
的 可 行 性 。 

总 体 说 来 ,80 年 代 关 于 双 波 长 方法 的 研究 主要 在 理论 阶段 ,包括 敏感 波段 的 选择 , Cy, 的 取 值 ,空气 折射 
率 结构 参数 与 温度 湿度 结构 参数 的 关系 等 方面 。 在 应 用 过 程 中 ,荷兰 Wageningen 大 学 利用 双 波 长 方法 在 不 
同 地 形 条 件 下 测算 区 域 平 均 蒸 散 ,修正 了 部 分 计算 方法 ,推动 了 微波 闪烁 仪 的 发 展 ;荷兰 Eindhoven 大 学 设计 
了 波长 为 11mm ,频率 为 27GHz 的 微波 闪烁 仪 , 常 与 LAS 联合 应 用 在 多 个 陆 面 实验 中 ,瑞士 Bern 大 学 ,对 双 波 
长 测算 方法 和 过 程 中 的 关键 参数 进行 研究 ,提高 了 测算 精度 5 ;德国 Radiometer Physics GmbH 公司 率先 研发 
了 波长 为 1.86mm ,频率 为 160GHz 的 商业 化 微波 闪烁 仪 ,为 双 波 长 的 发 展 提供 了 硬件 支持 。 


3 光 闪 烁 方法 测量 区 域 蒸 散 的 不 确定 性 


3.4 仪器 自身 限制 所 引起 的 不 确定 性 

光 闪 烁 方法 利用 测 得 的 C? 结合 MOST 进而 实现 对 区 域 蒸 散 的 估算 。 因 此 光 闪 烁 方法 的 仪器 ,如 LAS, 
MWS 等 ,会 因 其 自身 的 局 限 性 ,如 孔径 尺寸 .响应 时 间 等 ,在 大 气 平稳 状态 或 庙 流 很 强 的 状态 都 会 造成 波动 信 
号 的 偏差 ,其 中 最 显著 的 问题 就 是 信号 饱和 。Kohsiek 4&7 对 LAS 研究 发 现 ; 发 生 饱 和 现象 时 会 使 观测 结 
偏 低 约 20% 。 也 有 研究 表明 LAS 观测 数据 的 误差 会 随 信 号 强度 的 降低 而 增 大 "目前 针对 饱和 信号 的 处 
理 方法 已 有 不 同方 法 进行 研究 ,如 将 0.193552 A D? (L 为 路 径 长 度 以 为 闪烁 仪 所 用 波长 ,D 为 闪烁 仪 孔 径 
长 度 ; 下 同 ) 作 为 信号 阔 值 ,也 有 人 研究 表明 信号 羡 值 强度 为 0.074L °A O D? a, 0.099 L ?? A'2 p?? ,可 采用 饱和 
修正 因子 的 方法 或 几何 路 径 长 度 方法 对 饱和 信号 进行 修正 ,但 这 些 人 研究 均 存 在 不 足 , 所 得 结果 也 有 差 
异 *”]。 水 汽 也 会 对 光 信 号 产生 吸收 效应 ,由 于 此 效应 的 存在 , 光 闪 烁 方法 的 观测 结果 一 般 存 在 10% 左 右 的 
误差 ,极端 条 件 甚至 能 达到 80% ,但 是 由 于 水 汽 吸 收 引起 的 误差 可 通过 适当 的 滤波 方式 来 降低 六 。 

光 闪 烁 方法 的 仪器 也 存在 系统 误差 。 以 LAS XH, m= Kipp&Zonen 公司 同型 号 闪烁 仪 间 的 观测 误差 可 
达 2196 ,德国 Scintec 公司 同型 号 闪烁 仪 观测 误差 为 3% 左 右 !550 。 刘 绍 民 等 259 在 密云 站 的 比较 也 得 出 LAS 
系统 间 存 在 误差 ,并 且 试 验 中 发 现 Kipp&Zonen 的 LAS 结果 比 Scintec 的 BLS 型 号 闪烁 仪 观测 的 结果 存在 
17% 左 右 的 误差 。 

双 波 长 方法 采用 两 套 不 同 波长 的 仪器 协同 工作 ,为 获得 较为 可 靠 的 信号 ,通常 使 用 高 频 信号 观测 然后 进 
行 滤 波 处 理 。 然 而 不 同 研 究 者 ,所 使 用 的 频率 ,采用 的 波长 各 不 相同 ,如 表 1 所 示 。 


R1 双 波 长 方法 中 微波 选择 


Table 1 Selection of micro-wave in two-wave-length method 


作者 年份 Author, Year 波长 Wave-length/mm 频率 Frequency/ GHz 实验 地 点 Set-up site 
Kohsiek 等 (1983》 [40] 10 30 Eindhoven , fnj = 
Hill 4&( 1988)! 1.7 173 Flatville , 美国 
Green 等 (2001) [57 11 27 Kaitaia ,新 西 兰 
Lüdi 等 (2005) !'*! 3.2 94 Lindenberg ,德国 
Meijninger 等 (2006) [53] 11 21 Flevoland ,荷兰 
Evans( 2009) [53] 3.2 94 Sheepdrove , 英国 
Ward 等 (2015) 14 3.2 94 Swindon , 英国 
Mei-Sun Yee 等 (2015) P% 3.2 26,38 Yanco ,澳大利亚 


WRR KAE PEXEZEHIUAS REFUS b E RETE, JU Andreas! 在 研究 中 指出 ,波长 为 3.33mm 
时 , 波 文 比 在 0.1 一 10 范围 内 ,估算 的 蒸 散会 有 20% 的 误差 , 当 波 文 比 在 0.01 一 0.5 时 ,结果 会 有 10% 的 偏差 。 
3.0 ”环境 因素 引起 的 不 确定 性 

蒸 散 是 水 汽 传输 过 程 的 反映 ,准确 获取 蒸 散 量 主 要 取决 于 对 空气 中 热量 和 水 汽 含量 的 精确 测定 。 波 文 比 
为 显 热 通 量 与 潜 热 通 量 的 比率 ,主要 与 空气 温 湿度 有 关 , 可 以 用 来 表征 周围 环境 的 干 湿 状 况 。 单 波长 方法 的 
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人 研究 中 ,使 用 波 文 比 来 表征 水 汽 含量 ,在 不 改变 计算 中 的 其 他 参数 时 , 波 文 比 系数 由 1 变 成 0.3 时 可 造成 15% 
的 误差 :304 。 研 究 也 表明 , 波 文 比 在 0.1—2 之 间 变 化 时 结果 的 标准 差 在 0.4 左右 变化 ,对 结果 的 不 确定 性 影响 
可 达 5096 ,并 认为 当 波 文 比 大 于 0.32 时 ,对 结果 至 少 产 生 10% 的 影响 ,由 此 建议 用 饱和 水 汽 压 来 表征 空气 湿 
度 ' 5 Lüdi 等 (5 在 研究 中 用 绝对 湿度 表示 空气 湿度 ,结果 表明 当 波 文 比 系数 小 于 -2 时 ,相关 性 低 至 -0.6， 
而 当 波 文 比 系数 大 于 2 时 ,相关 性 为 0.8。 双 波长 方法 中 研究 人 员 多 倾向 使 用 绝对 湿度 ,而 Ward 等 5 认为 绝 
对 湿度 受 温度 影响 较 大 ,是 变量 ; 比 湿 独立 于 温度 ,相对 保守 ,并 在 Swindon 的 城郊 中 用 绝对 湿度 研究 ,结果 与 
涡 动 相关 方法 仅 相差 8% 5 。 在 光 闪 烁 方法 中 ,如 何 真实 地 表示 空气 湿度 状况 还 有 待 进行 更 深入 的 研究 。 

满足 相似 理论 (MOST) 是 兴办 烁 方法 测量 区 域 莹 散 的 前 提 。 不 同 的 植被 覆盖 面积 .土地 利用 方式 和 地 形 
起 伏 等 因素 造成 了 实际 观测 区 域 的 非 均 匀 性 , 光 径 高 度 随地 形 起 伏 不 断 变化 。 为 解决 MOST 在 不 均匀 地 表 的 
适用 性 问题 ,提出 了 “ 摊 混 高 度 ” 概 念 ,并 认为 在 此 高 度 条 件 下 所 有 通 量 能 够 充分 混合 '“1 ;而 Hartogensis 55:77 
人 研究 认为 ,在 不 稳定 状态 下 以 及 MOST 有 效 时 ,光路 径 上 的 通 量 值 主要 受 路 径 权 重 函 数 的 影响 ,提出 “有 效 高 
度 ” 的 概念 ,并 认为 当 路 径 距 离 大 于 5km 时 ,由 于 地 球 曲面 的 影响 会 产生 0.5m 的 偏差 ;Evans 58 研究 “ 双 波 
长 "系统 的 有 效 高 度 , 认 为 在 理想 的 下 垫 面条 件 下 有 效 高 度 引 起 的 误差 不 足 1% ,但 在 地 形 起 伏 较 大 的 情况 下 
有 5% 一 10% 的 误差 , 且 路 径 权 重 函数 对 路 径 终端 的 影响 也 较为 显著 。Ward 等 (中 认为 不 同 的 有 效 高 度 对 结 
果 产 生 6% 或 3% 的 偏差 。 地 形 起 伏 导 致 光 径 高 度 发 生变 化 ,MOST 的 表达 形式 因此 受到 影响 ,虽然 大 量 实验 
得 到 的 观测 资料 计算 分 析出 的 函数 表达 形式 基本 相似 (公式 (18) (19) ) ;但 是 函数 表达 式 的 具体 参数 仍然 相 
差 很 大 。 由 于 满 流 运动 及 其 影响 机 制 的 复杂 性 ,在 不 同 条 件 下 的 观测 实验 仍然 无 法 得 到 一 致 的 函数 表达 式 的 
具体 参数 ,如 表 2 所 示 。 


X2 不 同 研究 中 的 MOST 函数 系数 取 值 
Table 2 Coefficients of MOST function in different research 
MOST 函数 系数 Coefficients of MOST function 


系数 Coefficients 不 稳定 状态 Unstable 稳定 状态 Stable 

" " N 4.2 "n X2 xm x 
Wyngaard 等 ( 1973) 158] 4.9 7 — 一 4.9 2.75 一 一 
Andreas (1988) [59] 4.9 6.1 一 一 4.9 2.2 
Hill 等 ( 1992 ) [60] 8.1 15 一 一 一 一 一 一 
Hartogensis 等 (2005) [611 4.9 9 一 一 4.7 1.6 一 一 
Li & (2012) [€] 6.7 14.9 3.5 4.5 4.5 1.3 3.5 2.4 
XE (2013) [83 4.91 6.51 — — 6.43 2.2 — — 
Ward 45 (2015) 14 4.84 11.6 4.84 11.6 4.9 2:5 4.9 2.2 


MOST: 相 似 理 论 , Monin-Obukhov Similarity Theory ; xz 1 XT,2 54,1 34,3 : MOST PRICES 8C, coefficients of MOST function 


不 同 的 普 适 函数 关系 式 对 蒸 散 结果 影响 可 达 1096—1596 ,并 认为 在 不 稳定 状态 下 MOST 的 适用 性 比 在 近 
中 性 或 稳定 状态 时 好 , 汕 流 越 强 ,适用 性 越 好 !2 。 由 于 MOST 要 求 观测 试验 尽 可 能 处 于 均一 稳定 的 观测 条 
件 下 进行 ,但 当 热 力 庙 流 的 作用 减 小 , 主要 是 机 械 剪 切 力 、 重 力 及 平流 的 作用 ,MOST 适用 条 件 在 此 种 稳定 度 
条 件 下 不 再 满足 ,这 也 导致 了 不 同 普 适 函数 计算 结果 差异 很 大 。 闪 烁 仪 测量 的 是 大 气 油 流 强度 ,并 不 能 判断 
水 热 通 量 的 输送 方向 。 大 气 稳定 度 的 判断 及 MOST. 函数 的 选取 直接 影响 到 闪烁 仪 观测 的 数据 质量 ,因此 许多 
研究 人 员 只 分 析 白 天 观测 数据 ,很 少 考虑 夜间 观测 结果 。 目 前 对 于 普 适 函数 的 选择 尚 无 统一 标准 ,尤其 利用 
交 闪 烁 方法 计算 稳定 条 件 下 的 通 量 还 存在 较 多 问题 。 
3.3 ”计算 过 程 的 不 确定 性 

利用 单 波长 方法 ,通过 能 量 平衡 的 方法 计算 蒸 散会 造成 低估 ,即使 利用 波 文 比 修正 方法 进行 湿度 订正 后 ， 
结果 仍 有 1696 85/2: 77! ,其 原因 通常 认为 是 用 能 量 平衡 余 项 法 时 参与 计算 的 净 辐 射 (R,) 和 土壤 热 通 量 (C,) 
的 空间 代表 性 不 足以 代表 测量 区 域 的 均值 。 在 利用 能 量 平衡 余 项 法 时 , 净 辐 射 (R,) 和 土壤 热 通 量 ( G.,) 是 非 
常 重要 的 参数 ,在 非 均 匀 下 热 面 时 的 空间 差异 性 需要 仔细 测算 。 在 复杂 下 热 面 条 件 时 ,粗糙 度 长 度 、 零 平面 位 
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移 以 及 光 径 有 效 高 度 等 重要 参数 会 随 着 下 垫 面 的 实际 情况 而 发 生变 化 ,存在 不 确定 性 。 

应 用 双 波 长 方法 中 的 “Hil 双 波 长 "方法 时 ,大 多 数学 者 假设 Rn =+1。 实 际 研究 表明 , 当 Rj, =+1 时 ,其 
结果 不 能 很 好 地 反映 半 干 旱地 区 的 藻 散 状况 ,研究 发 现 ,不 稳定 状态 时 ,Rnw 取 +0.7 至 +0.95 范围 更 适 
合生 1。 这 种 假设 对 于 非 理 想 条 件 并 非 总 是 成 立 ,因此 Ludi 等 提出 实时 获得 Rj 的 交互 方法 ,并 认为 当 环境 
波 文 比 系数 范围 变 小 时 ,Rj 的 范围 也 会 缩小 ,并 认为 Rj 取 值 会 对 燕 散 结果 产生 20% 一 50% 的 偏差 。“Liidi 
双 波 长 方法 ”获得 实时 Rn 的 值 ,但 该 方法 中 所 用 的 公式 (8) 中 的 系数 4 要 依据 具体 安装 方案 来 定 ,仪器 的 间 
距 以 及 光路 径 长 度 的 测量 误差 会 引起 该 系数 的 变化 ,同时 该 方法 需要 近 红 外 路 径 与 微波 路 径 的 交点 位 于 路 径 
中 点 位 置 ,而 实际 测量 中 很 难 满足 此 要 求 。 

3.4 通 量 源 区 评价 的 不 确定 性 

目前 所 使 用 的 通 量 源 区 模型 对 平坦 均匀 下 垫 面条 件 适 用 性 强 , 在 复杂 的 下 垫 面条 件 下 则 有 许多 限制 ,如 
有 效 高 度 的 限制 ,粗糙 度 的 变化 .大气 稳定 度 的 变化 等 '““; 。2002 年 Meijninger “首次 将 源 区 模型 应 用 至 光 
闪烁 方法 测量 结果 的 分 析 ,并 认为 由 于 源 区 的 不 同 ,观测 结果 存在 8% 左 右 的 误差 全 区域 尺度 范围 内 ,由 于 通 
量 源 区 的 变化 导致 蔡 散 估算 偏差 会 高 达 2696 , Ward. 等 ' 和 1 也 认为 由 于 源 区 的 改变 而 造成 满 流通 量 的 低估 
也 是 引起 测量 结果 不 准确 的 原因 。 青 藏 高 原 不 同 空间 尺度 的 研究 表明 ,风向 和 通 量 源 区 的 不 同 会 导致 结果 的 
差异 ,并 说 明 通 量 源 区 不 重合 时 , 涡 动 相关 方法 与 LAS 观测 结果 差异 显著] ,对 蒸 散 的 影响 最 大 可 达 23% 左 
右 。Evans 和 De Bruin' 引 认为 使 用 源 区 函数 对 结果 进行 评估 时 必须 考虑 源 区 函数 模型 的 适用 条 件 并 指出 
Schuepp 模型 只 适合 在 平坦 下 垫 面 ,复杂 起 伏地 形 条 件 下 风向 会 引起 较 大 的 误差 , 源 区 函数 经 过 修正 后 得 到 
的 结果 更 好 :0 。 应 用 源 区 模型 应 考虑 观测 高 度 < 风 速 风向 、 下 扑面 的 粗糙 程度 以 及 大 气 稳定 度 等 因素 影 
响 ,不 同 气象 条 件 下 的 源 区 分 布 不 同 ,并 且 源 区 内 的 通 量 贡 献 值 各 不 相同 ,在 非 均 匀 下 热 面条 件 下 还 应 依据 源 
区 函数 对 最 终结 果 进 行 修 正 。 光 闪烁 方法 测量 区 域 尺度 通 量 时 ,由 于 下 垫 面 植被 类 型 不 同 , 各 斑 块 地 段 对 通 
量 贡献 也 不 相同 ,在 评价 过 程 中 需 进 一 步 探讨 。 此 外 ; 双 波 长 方法 中 由 于 两 套 系统 的 分 离 ,在 一 定 程度 上 会 造 
成 二 者 源 区 的 不 同 ,目前 关于 双 波 长 方法 的 源 区 仍 参考 单 波长 方法 ,这 也 增加 了 结果 的 不 确定 性 。 


4 结论 与 展望 


应 用 光 闪 烁 方法 的 实验 结果 表明 5; 光 闪 烁 方法 可 以 获得 区 域 水 平 上 的 平均 水 热 通 量 , 由 于 平均 时 间 比 涡 
动 相 关 方 法 短 ,因而 其 统计 不 确定 性 小 z 下 垫 面 均匀 时 , 光 闪 烁 法 的 测量 结果 与 涡 动 相关 方法 测量 结果 十 分 
一 致 ;而 在 非 均匀 下 垫 面条 件 下 ,充分 考虑 风向 与 通 量 源 区 的 影响 后 , 光 闪 烁 方法 也 可 以 取得 很 好 的 结果 ,万 
其 是 与 调动 相关 方法 联合 使 用 不 仅 能 获得 更 好 的 区 域 平 均 蒸 散 结果 ,而 且 也 能 精确 评估 各 源 区 的 蒸 散 水 平 。 
综合 看 来 , 兴 内 烁 方法 具有 如 下 优点 : 

(1) 快 速 测算 出 公里 尺度 路 径 上 的 平均 水 热 通 量 ,在 复杂 下 垫 面条 件 下 更 具 优 势 ; 

(2) 是 目前 卫星 遥感 模型 反 演 验证 的 最 佳 地 面 实测 方法 ; 

(3) 具 有 较 高 的 时 间 分 辩 率 ,空间 代表 性 强 ; 

(4) 具 有 在 复杂 天 气 条件 下 全 天 候 连续 观测 的 能 

就 光 闪 烁 方法 而 言 , 单 波长 方法 发 展 更 为 成 熟 ,已 经 完全 实现 了 商业 化 。 单 波长 方法 测算 蒸 散 需 额 外 配 
备 净 辐射 .土壤 热 通 量 等 能 量 传 感 器 ,通过 能 量 平衡 公式 计算 得 出 获 散 ”人 ,但 其 中 过 程 变量 如 净 辆 射 E 
壤 热 通 量 等 与 光 闪 烁 方法 的 测量 路 径 ( 公里 尺度 ) 存 在 空间 差异 性 。 单 波长 方法 选用 的 近 红外 波段 会 受 
到 湿度 的 影响 ,计算 时 通常 参考 Wesely 引入 的 波 文 比 修正 系数 ‘3 ,并 没有 考虑 实时 的 湿度 变化 ,通过 能 量 平 
衡 计 算 蒸 散 的 精确 性 尚 需 提 高 ?1 。 双 波长 方法 同时 测量 观测 区 域 的 温 湿度 变化 ,一定 程度 上 消除 了 单 
波长 方法 计算 蒸 散 过 程 中 存在 的 空间 差异 性 以 及 湿度 变化 的 影响 。 虽 然 双 波长 方法 的 发 展 起 步 较 晚 ,研究 报 
道中 所 用 的 双 波 长 仪器 大 多 为 实验 室 自行 研发 ,商业 化 仪器 较为 少见 ,价格 较为 昂贵 ,普及 程度 不 及 单 波 长 方 
法 的 仪器 ,但 是 双 波 长 方法 简化 了 蒜 散 测算 的 过 程 ,实现 了 对 蒜 散 的 直接 测量 ,在 未 来 高 精度 的 区 域 通 量 研究 
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中 具有 更 大 的 优势 。 


总 之 , 光 内 烁 方法 在 短 短 几 十 年 间 得 到 了 迅速 发 展 ,并 成 为 区 域 落 散人 研究 领域 的 重要 方法 ,特别 是 双 波 长 
方法 在 非 均 匀 地 表 燕 散 观 测 研 究 中 表现 出 了 巨大 的 发 展 潜力 。 光 闪烁 法 不 仅 实现 了 区 域 通 量 的 连续 观测 , 同 
时 也 验证 了 遥感 反 演 结果 ,改善 了 遥感 模型 参数 ,为 区 域 尺度 上 的 水 资源 管理 做 出 了 重要 贡献 ,但 光 闪 烁 方法 
还 存在 如 信号 饱和 、 相 似 函 数 选择 等 不 确定 性 。 如 何在 最 大 程度 上 减少 不 确定 性 ,作者 认为 可 以 从 下 述 几 点 
进行 进一步 的 研究 。 

(1) 光 闪烁 方法 在 理论 方面 的 研究 。 主 要 包括 光 闪 烁 方法 在 测算 由 于 频率 损耗 等 原因 引起 的 灌流 谱 修 
正 , 饱 和 效应 修正 ,水 汽 吸 收效 应 修正 ,结合 清流 谱 分析 与 涡 动 相关 方法 观测 构造 合适 的 相似 函数 。 

(2) 计算 过 程 中 关键 参数 的 确定 以 及 源 区 分 析 。 主 要 包括 复杂 下 垫 面条 件 下 粗糙 度 SEEDS. \ 有 效 高 
度 等 参数 的 确定 ,复杂 地 形 条 件 下 的 源 区 分 析 等 。 

(3) 其 他 方面 的 研究 。 如 光路 径 上 存在 的 气 洲 胶 吸 收 影响 、 双 波长 方法 中 仪器 分 离 造 成 的 路 径 权 重 函 数 
分 离 . 光 闪 烁 方法 测算 结果 的 验证 等 。 
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